COMPETITIVE EXCLUSION FOR CHEMOSTAT EQUATIONS WITH 

VARIABLE YIELDS 

TEWFIK SARI 

Abstract. In this paper, we study the global dynamics of a chemostat model with a single 
nutrient and several competing species. Growth rates are not required to be proportional to food 
uptakes. The model was studied by Fiedler and Hsu [J. Math. Biol. (2009)59:233-253]. These au- 
thors prove the nonexistence of periodic orbits, by means of a multi-dimensional Bendixon-Dulac 
criterion. Our approach is based on the construction of Lyapunov functions. The Lyapunov func- 
tions extend those used by Hsu [SIAM J. Appl. Math. (1978) 34:760-763] and by Wolkowicz and 
Lu [SIAM J. Appl. Math. (1997) 57:1019-1043] in the case when growth rates are proportional 
to food uptakes. 
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1. Introduction and main result 



The aim of this paper is to prove the Competitive Exclusion Principle (CEP) for the following 
competition for a single resource model 

N 

(1) S' = D(S°-S)-Y, Pi (S) Xi 

i—\ 

x'i = fi(S)x iy i = l---N, 

where S(t) and Xi(t), i = 1 • • • N, denote respectively the nutrient concentration and the concen- 
tration of the ith competing species at time t. The input concentration 5° and the removal rate D 
are assumed to be constant. The uptake rate Pi(S) satisfies 

O ' ( 2 ) Pi(fy = and Pi(S) > for S > 0. 



The growth rate fi(S) satisfies 

(3) MS) < for < S < A, and /.(A,) = 0. 

The smallest positive value of the concentration substrate S — \ given in §3§ , where the growth 



rate fi(S) of Xi is 0, is called the break-even concentration for the ith species. This model was 



considered by Fiedler and Hsu [7] as an extension of the classical chemostat model 

N 

S' = D[S°-S}-J2- 






(4) u ^ ^ ^ Y t 

i— l 

x[ = [qi{S) — Di]xi, i = l---N, 
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2 TEWFIK SARI 

where Yi is the conversion constant, or yield, for the ith species, and the growth rate qi(S) satisfies 

(5) g;(0) = and qi (S) > for S > 0. 

Thus, the break-even concentration 5* = A; satisfies the condition q%(S) = Df, it is the smallest 
value, where the growth qi(S) of x% is balanced by the removal rate Di. The chemostat occupies 
a central place in mathematical ecology, see the monograph of Smith and Waltman [23]. It is a 
model of the dynamics of microbial competition. Basically, the chemostat consists of a nutrient 
input, pumped at a constant rate into a well- mixed culture vessel. The culture vessel contains 
the microorganisms that are growing and competing for the nutrient. Volume is kept constant by 
pumping the mixed contents out at the same rate. If the growth functions fi(S) and the uptake 
functions Pi{S) are of the form 

(6) MS) = q t (S) - D u Pi (S) = qi(S)/Yi, 

then the model (TTJ) of Fiedler and Hsu [7J reduces to the simple chemostat model ((H). However, 
the particular form ^ of the growth function fi(S) is not assumed in [7] and hence, the model (fTJ) 
recovers the more general case of variable yields model 

N 

(7) S' = D[S°-S\-J2Pi(S)xi, 

x\ = [ ft (5) - D t ]xu i=l---N, 

where uptake rate p, (5) satisfies @ and growth rate qi(S) satisfies §5§. The functions yi(S), defined 
by 

(8) y i (S) = ^r, i = l-N, 

Pi{S) 

are the growth yields. The model (J7J) was considered by Arino, Pilyugin and Wolkowicz [2J as a 
generalization of the constant yields case (UJ), for which the yields © satisfy yi(S) = Yi. It was 
further investigated by Sari and Mazenc [21] [22] • Without loss of generality (see Section [2]), we 
assume that D = 1 and S a = 1 in ([!]). The system becomes 



JV 

a' — 1 a 

(9) 



S' = l-S-J2Pi( S ) x i 



x'i = fi(S) Xi , % ' i = l---N. 

Coexistence of the N species is a fundamental question on the model ^ of competition for a 
single resource. Looking for coexistence at positive equilibria we have to solve equations fi(S) = 
simultaneously for all i = 1 • • • N. In general, for N > 2, these equations cannot be solved for the 
same value of S. Thus, generically, §§§ can have the following equilibria: the washout equilibrium 

(10) £ = (1,0,---,0), 

where all species go extinct, and equilibria Ei, i — 1 • ■ • N, where all components of Ei vanish, 
except for the first and the (i + l)th, which are 

1 — S 1 * 
S = S* i 



Ms*) 1 

where S* S]0,1[ satisfies fi(S*) = 0. Hence, at any equilibrium point Ei, all but one species go 



extinct. 
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Since /(Aj) = 0, the break-even concentration S* — Xi gives rise to an equilibrium point Ei 
for the system, if and only if A,; < 1. A well-known open-problem in the theory of the chemostat 
is to prove the global asymptotic stability of the equilibrium point Ei with the lowest break-even 
concentration. If this equilibrium is globally asymptotically stable (GAS), then the CEP holds: 
only one species survives, namely the species which makes optimal use of the resource. The reader 
is referred to [24J, for complements and details on the CEP, and to [TJd] for recent results and a 
discussion on competitive exclusion. Most of the results on the CEP for (0| and (|7|) have been based 
on Lyapunov functions (3j [9j [14j [2TJ [22J, [26l [27] . For a survey of constructing Lyapunov functions 
in the chemostat, the reader is referred to [10] . We simply recall here that Hsu |9J proved the CEP 
for the Monod case of ((4]), when the growth functions are 

(ID *< s > = ^ 

and Wolkowicz and Lu [2B] extended the result of [5] to more general growth functions. 

Instead of a Lyapunov function approach, Fiedler and Hsu [7] applied a multi-dimensional 
Bendixon-Dulac criterion to exclude periodic solutions. Under some technical conditions on the 
functions fi and pi they proved that ((HJ) does not possess positive non-stationary periodic orbits. In 
our previous works [21, 22J, we showed that both Lyapunov functions used by Hsu [9j and Wolkow- 
icz and Lu [55] can be extended to the variable yields case model ([7]) . The aim of this paper is to 
show that these Lyapunov functions can also be used to obtain the CEP for ©. 

A necessary condition to avoid washout of the species, and global convergence towards the 
washout equilibrium E defined by (flu) , is that Aj < 1 for at least one species. Assume that the 
species are labeled so that < Ai < 1. Then 

(12) El = (X x ,x$,0,- ■■,()), 

where x\ — x\ = Pi(Ai) is an equilibrium. Here 

(13) P ^ S )= X —uk- 

Pi(S) 

Using linearization of ^j about E* one proves that: 

Lemma 1.1. The equilibrium Itlty) is locally exponentially stable if and only if /{(Ai) > and 
P{(Xi)<0. 

We consider the global asymptotic stability of Ef . Our main result is 

Theorem 1.2. Assume that Ai < 1 and for all < S < 1, 

(14) (S-X 1 )f 1 (S)>0, forS^X,, 

(15) (S - X 1 )(P 1 (S) - i^(Ai)) < 0, for S ± X x . 

Assume that there exist constants ai > for each i > 2 satisfying A; < 1, such that for all < S < 1, 

(16) f 1 (S) Pl (S)>a l f l (S)(l-S). 

Then the equilibrium E£ is GAS for (0) with respect to the interior of the positive cone. 

Notice that the following property holds. 
Lemma 1.3. The conditions X± < 1 and Mb}) imply that Xi < A^ for all i > 2. 
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Proof Assume that there exists i > 2 such that A^ < Ai- Then, there exists S < Ai such that 
fi(S) > 0. Hence, since S < Ai < 1, /i(5)(l - S) > 0. On the other hand, using ©, /i(S) < 0. 
Hence, the inequality (p~6|) is violated. 

This lemma shows that the winning species x\ of Theorem 11.21 has the lowest break-even concen- 
tration, in accordance with the CEP for models of competition for a single resource |24| . 

The paper is organized as follows. In Section [2] we give some preliminary lemmas. In Section 
[3] we show how the Lyapunov function of Wolkowicz and Lu [26] can be extended to (j9]) and used 
to obtain Theorem 11.21 We show in this section that the result of [2 lj for 0, which extends the 
result of [26J for (01, is a corollary of Theorem 11.21 We give also graphical interpretations of the 
conditions (|14l) . (Tl5|) and (p~6|) in Theorem 1 1.2 1 In Section [4j we show how the Lyapunov function of 
Hsu [9] can be extended to (j9|) and used to obtain Theorem 14. 1[ which is another global asymptotic 
stability result of E*{ for ([9]) . Theorem 14.11 can be obtained also as a corollary of Theorem 11.21 (see 
Proposition 14. 2p . We show in this section that the result of [5] for (H|) with Monod functions (|lll) . 
and the result of [35] for ([7]) are corollaries of Theorem l4.ll In Section [5] we discuss the single species 
case N = 1. In Section [6] we apply our result to the model with Monod growth functions (fTTj) and 
linear yields. In Section [7] we discuss some of the CEP results based on Lyapunov functions and we 
compare Theorem 11.21 with the results of [7] based on a Bendixon-Dulac approach. 

2. Preliminary results 

Let us prove first that we can we assume that D = 1 and S° = 1 in ([I]). Indeed, under the 
change of variables 

-Q ^ T 7-,. - f-Q\ Pi(b S) -j .-pp. ji(b S) 

S =gQ, t = Dt, Pi (S) = SQD , fi{S) = — ^ — , 
equations (JTJ) take the form 



dS 1 dS 



dt S°D dt 

dxt 1 dxi 



N 

1 - S -^2pi(S)xi 



fi(S)xi, i = l---N. 



dt D dt 

Dropping the bars, one obtains ©. Recall that /,(0) < 0, so that the concentration of the species 
xi is decreasing when the concentration of nutrient is too small. The smallest positive zero S = \{ 
of fi is the break-even concentration of the ith species Xi. We adopt the convention A^ = oo if 
fi(S) < for all S > 0. We need the following lemmas. 

Lemma 2.1. The non-negative cone is invariant under the flow of |P|) and all solutions are defined 
and remain bounded for all t > 0. 

This lemma is simply Theorem 4.1 in [2]. 

Lemma 2.2. If for some species Xi, the inequality (S — Xi)fi(S) > is satisfied for all < S < 1, 
S J^ Xi, then S(t) < 1 for all sufficiently large t and all initial condition . 

This lemma can be obtained using arguments similar to that given in the proofs of Lemma 2.9 
in [5] and Lemma 2.1 in [26J. 

Lemma 2.3. For all solutions of (0), if Xi > 1 then Xi(t) — > as t -^ oo. 

This lemma can be obtained using arguments similar to that given in the proofs of Lemma 4.2 
in [5] and Lemma 2.2 in |26| . 
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3. Extension of the Lyapunov function of Wolkowicz and Lu 
The Lyapunov function used by Wolkowicz and Lu [26] in the constant yields case Q is 

(1?) VWL = "A- L -^^ d ° + Y X ] xl —^ + g Y*- 

with suitable constant a > 0. Using the notations in ©, and since S° was rescaled to 1, the 
numerator in the first integral of (fl7f is simply equal to ^_J . Multiplying (|17p by the constant 
TTx - = P"> §i ves the following function 

(18) V = I" f^da + -1 f X ^d£ + f; a lXu 

where a, are constants to be determined. This is a Lyapunov function for (JSJ) which permits to 
prove Theorem 11.21 as shown below. 

Proof of Theorem \1.2\ From Lemmas 12.21 and 12.31 it follows that there is no loss of generality to 
assume that Aj < 1 for i = 1 • • • N and to restrict our attention to < S < 1. Consider the function 

V — V(S,xi, ■ ■ ■ ,Xn) given by (|18|) where on are positive constants satisfying (|16|) . The function 

V is continuously differentiable for < S < 1 and Xi > and positive except at point E*. The 
derivative of V along the trajectories of ^ is 



h(S) a i , - T i 



v 



1 — b x\x\ *-^ 



OiiX 



l^i- 



=2 



Since x \ = Pi(Ai) and using ©, V is written 



y,_h{S) 



N 

l- s - y~]pi(s)xi 



=i 



l-S 

Y=g-0- —") """ Pi(Ai) 



1 



V 



, ■[x 1 -P 1 (A 1 )]/ 1 (5) + ^a i / i (5K 



The terms vie (1 — 5) and — p h ■. Pi(Ai)/i(6 f ) are canceled. Hence, using (|T^|) . 



^'=3:1/1 (5) 



1 1 



Pi(Ai) Pi (5) 



A a i / i (S)(l-S)-/i(5)p < (S) 
V* 1-^ • 



Using p^|) and (fT3|) , the first term of the above sum is non-positive for < S < 1 and equals if 
and only if S = Ai or x\ — 0. Using (fT6"|) , the second term is non-positive for < 5 < 1 and equals 
if and only if Xi = for i = 2 • • • N. Hence V < and V = if and only if Xi = for i = 1 • • • N 
or S — Ai and a;^ = for i = 2 ■ ■ ■ N. Using the Krasovskii-LaSalle extension theorem, the w-limit 
set of the trajectory is E^. 

Theorem 11.21 was previously obtained in |21j . in the particular case when the function fa has at 
most two positive zeros A^ and /ij, with Ai < /^ < +00, such that 

(19) MS) < if S i {Xi,Hi}, and f t (S) > if S E}\ u ^{, 

with the convention that \Xi = 00 if equation fi(S) = has only one solution and Ai = 00 if it has 
no solution. This class of functions corresponds to the case when fi(S) = qi(S) — Di and 

qi(S) < Di if S g [Xi,fj,i], and q^S) > D { if S e]Aj,^[. 
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It was often considered in the literature [H HH US] [57] . For this class of systems the main result in 
[21] is 

Corollary 3.1 (Theorem 2.1 in [21j). Assume that 

(20) Ai < A 2 < • • • < \n, and Ai < 1 < j«i, 

(21) (S - Xi)(Pi(S) - Pi(Ai)) < 0, for S^Xl 
There exist constants a > for each i > 2 satisfying Xi < 1, s-uc/i i/iai 

(22) max .gi(.S') < q < min gi(S), 

0<S<Ai \i<S<pi 

where gi(S) — 7a S hsmi-s) an< ^ P' 1 = mm (Mjj -0- ^ erl ^ e equilibrium E\ is GAS for @) wii/i 
respect to the interior of the positive cone. 

Proof Assume that ([20]). (|2"Tjl and ([22]) hold. Let us prove that ([14]), ([15]) and ([16]) hold. First, 
note that ([2"T]) is the same as ([T5]) . and condition Ai < 1 < /.ii in ([20]) is equivalent to ([T4]) . If 
Ai < S* < Ai then /i(S') < and /i(S') > so that ([TH]) is satisfied for any choice of ai > 0. 
Similarly if /i, < 1 and /x, < 5 < 1 then /^ (S 1 ) < and /i (S 1 ) > so that (|16p is satisfied for any 
choice of a. L > 0. On the other hand, if < S < Ai then fi(S) < and, using gi(S) < Ci in ([22]) . 

/ 1 (S) ft (5)>c i ^^/ < (S)(l-5). 
Finally, if A, < S < p i7 then fi(S) > and, using #;(£) > q in ([22]) . 

/ 1 (S) ft (5)>c i fi^/ i (S)(l-5). 
1 — Ai 

Thus ([T^]) is satisfied for o^ = c ^j 1 ^ 1 ^ . The result follows from Theorem II .21 

Condition ([T4"]) means that S 1 = Ai is the only zero of the growth function fi(S) for < S < 1. 
Condition (fT"5"]) means that S — Ai is the only zero of the function Pi(S) given by ([T5]) . for < S < 1. 
The technical condition ([T6]) is trivially satisfied in the single species N = 1. Following [22, 26] we 
give now a graphical interpretation of ([16]) . For each i > 2 such that A.; < 1, consider the function 



(23) 9l (S) = ^\ t f - 

fi{S)pi{S) 

The functions <7j is defined on (0, Ai) U (Ai, 1]. It tends to ±oo when S tends Ai. This function 
should not be confused with the function gi in Corollarv l3 . 1 1 which is simply a multiple of the inverse 
of gi. We use gi instead of its inverse, since the zeros of fi on [0, 1] are not known as for the class 
of functions fi considered in Corollary 13. II Since f\(S) < over [0,Ai) and fi(S) > over [Aj,l], 
the condition (|16p is equivalent to 

(24) min gJS) > — > max gAS). 

o<s<\i at \,<s<i 

Hence, the condition (TTOj) in Theorem 1 1.2 1 can be easily depicted graphically: plot simply the graph 
of gi(S) over [0, 1] and see if there is a gap between mino<s<A 1 9i(S) and max>. < s < i gi(S), see Fig. 
® 

It was shown in [ZT\ that the main result (Theorem 2.3) of [26J is a consequence of Corollary 13. II 
Hence, it is also a corollary of Theorem 11.21 
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4. Extension of the Lyapunov function of Hsu 

The Lyapunov function Vh used by Hsu [9] in the Monod case of (j4]), where the growth functions 
are of the form (fTTj) . is 



CiXi , 

with constants 

b 1 D 1 



(25) V H = I" G -^do + Cl r L-^ + V . 



AT, and Ai 



Fi Oj -A ax - D 



The function in the first integral of (|25)1 is simply given by „ Al = ci ; J , where 

oi + S Y\ bx + b 

Hence, multiplying (|25p by the constant 1/ci gives the following function 

(26) V = [ S ^da+ f^ + E«. 

where the constants a/ci in the last sum are simply denoted by a to avoid unnecessary new 
notations. Under some technical conditions, this function is a Lyapunov function for ((9j) and 
permits to obtain the global asymptotic stability of the equilibrium point El as stated in the 
following result. 

Theorem 4.1. Assume that Ai < 1 and for all < S < 1, 

(27) (S - Xx)fi(S) > 0, forS^Xx, 

(28) (S - Ai)(Pi(5) - Pi(Ai)) < 0, for S ^ \x. 

Assume that there exist constants c-i > for each i > 2 satisfying Ai < 1, such that for all < S < 1, 

(29) h(S) Pi (s) > «ms)pi(S). 

Then the equilibrium E[* is GAS for (0|) with respect to the interior of the positive cone. 



Proof From Lemmas 12.21 and 12.31 it follows that there is no loss of generality to assume that 
A-j < 1 for % = 1 • • • N and to restrict our attention to < S < 1 . Consider the function V = 
V(S,Xx, ■ ■ ■ ,Xn) given by (|26j) where a are positive constants satisfying (|29|) . The function V is 
continuously differentiable in the positive cone and positive except at point E*, where it is equal 
to 0. The derivative of V along the trajectories of (0) is 

v , _ h( S ) a' , x l - x *l ' , V r V 

Pi(S) xx f^ 



Since x\ = -Pi(Ai) and, using ([9]), V is written 



y, _ h(S) 

' Pi(S) 



N 



l-S~Y,Pi(S> 



N 

[xx-Pi(Xi)}h(S) + ^2cifi(S)s 

i=2 
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The terms — i L ) s l p 1 (S)xi and Xif±(S) are canceled. Therefore, using (fT3|). 



V' = fi(S)[P 1 (S)-P l (\ l )} + 



A Cifi(S) Pl (S) - fi(S) Pi (S) 

h Xi *(*) 



Using (f2"T)) and (j2"8)h the first term of the above sum is non-positive for < S < 1 and equals 
if and only if S = \\. Using (|29|) . the second term is non-positive for < S < 1 and equals if 
and only if Xi = for i = 2 • • • N. Hence V' < and V' = if and only if S — Ai and x% — for 
i = 2 ■ ■ ■ N. By the Krasovskii-LaSalle extension theorem, the w-limit set of the trajectory is E*. 

Notice that, as in Lemma [1.31 the condition (|29p implies that Ai < A.; for all i > 2. Indeed, 
if fi(S) > for some S < Ai, then /i(5) < 0, so that the inequality (|2"9"]l is violated. Thus the 
winning species x\ in Theorem 14.11 has the lowest break-even concentration. Actually Theorem 14. II 
is a consequence of Theorem 11.21 

Proposition 4.2. Theorem\4.1\ is a corollary of Theorem \l.Sl 



Proof Assume that flU]), (J2HI) and (gHJ) hold. Notice that (27]) is the same as CG]} and (JSHJ) is the 
same as (|15[) . Let us prove that (flo| holds. If fi(S) < and fi(S) > (which occurs if Ai < S < A; 
and may occur also for A, < S < 1), then (|16|) holds for any choice of ai > 0. If < 5 < Ai then, 
by (E5J), P^S) > Pi(Ai) and, since fi(S) < 0, ^§r > ^§j. Finally, if A, < 5 < 1 and f z {S) > 

then, by (J25J, Pi(S) < Pi(Ai), and hence, ^^ > jffif \ ■ Therefore, in both cases A 4 < 5 < 1 and 
0<5< Ai, 



M<?) > MS) 



Pl(S) Pi(Ai) 
Thus, using ([25]). 

/!(5) Pi (5) > aMS)MS) = c*^(! - 5 ) > C *;H§) C 1 - 5 )- 

Thus, (pi]) holds for a, = p^y- Hence, (J2SJ and fl2JjJ) imply fH]). This ends the proof. 

Theorem 14.11 recovers the classical case Monod case [9]. Indeed, consider the particular case of 
((3]), when the growth functions qi(S) are given by ((TT|) . System Q, with D = 1 and S° = 1, takes 
the form 



S' = 1 - S - £) 



TV „ 

Gin *J iCy 



(30) ~ bi ~' r S ' 



1 • ■ • N. 



M + s 

We consider the case where, for alH = 1 • • • N, ai > Dj. The break-even concentrations are 

(3D ** = ;r%- 

Corollary 4.3 (Theorem 3.3 in [9]). Assume that 

(32) Ai < A 2 < ■•• < X N , Ai < 1. 

TTien i/ie equilibrium E\ is GAS for Ii30\) with respect to the interior of the positive cone. 



COMPETITIVE EXCLUSION FOR CHEMOSTAT EQUATIONS WITH VARIABLE YIELDS 9 

Proof Assume that {3SJ holds. Let us prove that (J27J), (J2ZJ) and (J2ZJ) hold. Since f x (S) = 
qi(S) — D\ is increasing, the function fi(S) changes sign only at S = Ai and hence, ([27)) is satisfied. 
Since 

— and P 1 (S) = -Yi -^ 

the function Pi(S') changes sign only at S = Ai and hence ([55)) is satisfied. Condition 



J\(S)=y x (l-S)42L_ and P 1 '(5) = -F 1 — ^<0, 

I 

(ai -D^ff-fri-Di f aj5 (a. t - Dj)S - bjD., 1 a x S . 

h + S Y l b l + S >Ci h + S Yth+S' % ~ ' 

After simplification by ,, s)tb+s) » ^' s con dition is equivalent to 

(33) (a 1 -D 1 )^r(S-X 1 )>c i (a i -D i )^(S-X i ) i > 2, 

which is satisfied for C, = <' 'In'i V' • in dced, for this choice of the constants Cj, (|33p is simply 

S - Ai > 5 - Ai -«=>• Ai < A,, i>2, 

which is the same as (|32D . Thus (|2U1) is satisfied. The global asymptotic stability of E\ follows from 
Theorem 14.11 



Theorem 14.11 was previously obtained in [22], in the particular case when the function fi satisfies 
in]) . For this class of systems the main result in (52] is 



Corollary 4.4 (Theorem 2 in [22J). Assume that 

(34) Ai < A 2 < • • • < Ajv, and Ai < 1 < /zi, 

(35) (5 - A 1 )(P 1 (5) - Pi(Ar)) < 0, /or S ± X x . 
There exist constants ai > /or each i > 2 satisfying Ai < 1, suc/j i/iai 

(36) max ^(S 1 ) < c, < min Qi(S), 

0<S<Xi Xi<S<pi 

where gi(S) — /.(slisi an< ^ Pi = min(/ii, 1). Then the equilibrium E[* is GAS for ([I]) with respect 
to the interior of the positive cone. 

Proof First, note that (|35|) is the same as ([28)) . and condition Ai < 1 < /ii in (|34|) is equivalent 
to (|2"7), If 5 < Ai then / 4 (5) < and fi(S) > so that ([29]) is satisfied for any choice of a > 0. 
Similarly if \ii < 1 and /i^ < 5 < 1 then fi(S) < and fi(S) > so that (|29| is satisfied for any 
choice of a > 0. On the other hand, if < S < Ai then fi(S) < and, using gi(S) < c; in ([36]), 
fi(S) Pi (S) < Cifi(S) Pl (S). Finally, if A 4 < 5 < Pi , then / 4 (5) > and, using gi (S) > c, in flMD, 
fi(S)pi(S) < Cifi(S)pi(S). Thus {2JJ is satisfied. The result follows from TheoremO 

5. Single species 

In the case N = 1, @ takes the form 

, , S' = 1 -S-xp(S), 

[67) x' = /($>. 

Let 5 = A be the smallest positive value of S such that f(S) = and a;* = P(X) with P(5) 
defined by P(S) = ^-rM as in ([T3| . If A < 1, then _E* = (A, a;*) is a positive equilibrium. Assume 
that /'(A) > and P'(X) < 0, so that E* is locally asymptotically stable. We consider the global 
asymptotic stability of E* . 
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Corollary 5.1 (Theorem 2.11 in [2] or Lemma 2.3 in |TH])- Assume that A < 1 and for all < S < 1, 

(38) (S - X)f(S) > 0, for S ± A, 

(39) (S - X)(P(S) - P(X)) < 0, forS^X. 

Then the equilibrium E* is GAS for A37\) with respect to the interior of the positive cone. 



Proof Notice that ((38} is the same as flT4J) or ([27} and ((39} is the same as ([15} or ([28}. Since 
for TV = 1, condition ([TBI) in Theorem ll.2l or condition (|29p in Theorem 14. II is trivially satisfied, the 
result is a corollary of Theorem 11.21 or Theorem l4.ll 

Corollary 15.11 was obtained by Arino, Pilyugin and Wolkowicz (see [2], Theorem 2.11). Using the 
Lyapunov function 



(40) 



Vapw = 



1-A 
P(A) 



/to 



da 



i 



i 



-d£, 



these authors proved that if 

(4i) 1 *W0-x> 



has exactly one sign change for < S < 1 



p(X)(l-S) 

then E* is globally asymptotically stable. Notice that ([4"T} is equivalent to ([5T)} . In the single species 
case, our Lyapunov function ([T5} . used in the proof of Theorem 11.21 reduces (up to a constant) to 
the Lyapunov function Vapw considered in [2J . Corollary 15.11 was obtained also by Pilyugin and 
Waltman (see [18J, Lemma 2.3). Using the Lyapunov function 



(42) 



Vpw = 



m 



da 



? 



■«, 



IX PW " ' Jx* £ 

these authors proved that if 

(43) S = X is the only zero of R(S) = 1 - S - i*p(5) 

then S* is globally asymptotically stable. Notice that (l4"3"} is equivalent to (|39p . In the single 
species case, our Lyapunov function (|26p . used in the proof of Theorem l4.ll reduces to the Lyapunov 
function Vpw considered in |18j . 



x = P(S) 






A < Si or A > S 4 S 2 <X<S 3 Si < X < S 2 or S 3 < X < Si 

FIGURE 1. The graph of the function x — P{S) showing the values Si, S2, S3 and S4. 

Notice that the isoclines 5" = and x' = of ([37} are given by 

5" = ^=^ x = P(S), 
x' = <s=^ x = or S = X. 
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If the vertical line S = A intersects the curve x = P(S) on an increasing arc, then, from Lemma ll .11 
the intersection is an unstable equilibrium point E* . Using Poincare-Bendixon theory we can show 
that the system has at least a periodic orbit surrounding the equilibrium. Otherwise, if the vertical 
line S = X intersects the curve x — P(S) on an decreasing arc, then, from Lemma [L~T1 P* is locally 
asymptotically stable. The condition (|39|) has the following graphical interpretation: if E* is the 
only intersection of the isocline x = P(S) with the horizontal line x = x* then E* is GAS. For 
instance, in the situation depicted on Fig. [3J the function x = P(S) has two critical points S — S 2 
and S = S 3 . Let Si and S 2 defined by P(Si) = P(S 3 ) and P{S 4 ) = P{S 2 ) respectively. Then: 

Case 1. If < A < Si or S4 < A < 1 then E* is the only intersection of the isocline x = P(S) 
with the horizontal line x — x* . Thus, using Corollarv l5.11 the equilibrium E* is GAS. 

Case 2. If S2 < A < S3 then, using Lemma 11.11 the equilibrium E* is unstable. The system 
admits at least one limit cycle. 

Case 3. If Si < A < S2 or S3 < A < S4 then, using Lemma 11.11 the equilibrium E* is locally 
asymptotically stable. The horizontal line x = x* has three intersections with x = P{S). Since 
(f39l) does not hold, we cannot conclude if the equilibrium is GAS or not. 



We illustrate the third case by an example taken from [18J. Consider (|37l) with 

P(S) = 4fr and /( 5 ) = 9(S) - A where 9( 5 ) = 1^ and ^ 5 ) = l + c S 2 , 
y(b) b + S 

corresponding to Monod growth function and quadratic yield. For the parameter values given in 
the caption of Fig. [21 and D > D3 and close to D3, the equilibrium point E* is exponentially stable 
and it is surrounded by two limit cycles. Actually, the limit cycle which exists for all D2 < D < D3 
disappears for some critical D c > D3 through a subcritical Hopf bifurcation. For more details and 
explanations the reader is referred to |18j . 






S1S2 

FIGURE 2. If A < Si or A > S 4 , using Corollary 1 5. 11 the equilibrium point E* is 
GAS. If Si < A < S 2 or S 3 < A < S 4 the condition (|2H]) does not hold and Corollary 
15. II cannot be applied. Actually if D = 1 the system has two limit cycles. On the 
center of the figure an enlargement of the graph shows that the equilibrium point 
E* = (A, 2;*) lies on a decreasing branch of the graph of the function x — P(S). 
For the parameters values a = 2, b = 0.58 and c = 46: Si ~ 0.048, S 2 ~ 0.143, 
S3 ~ 0.579, S 4 ~ 0.855 and A = 0.58. Hence S3 < A < S 4 . 
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(44) 



The case N = 1 of a single species can also be investigated with the Bendixon-Dulac criterion. 
As shown in |7J, using the variable y — log (a;), ([37f is written 

S' = 1 - S - e*p(S), 

y' = f(S), 

with resulting divergence 

div = -1 - e y p'{S). 
If p'(S) < then the divergence is negative and no periodic solution can exist. Poincare-Bendixon 
theorem shows that convergence to the equilibrium E* ensues. Thus E* is GAS under the conditions 
p' > and (|38|) . This result is a consequence of Corollary 15.11 since the condition p'(S) implies 
P'(S) < for < S < 1 and hence, (J39J) holds. However, the condition ((39]) in Corollary 15. II can 
accept slightly negative p', since P'(S) < 0, for < S < 1, is equivalent to p'(S) > ~*_c i f° r 
0<S<1. 



6. MONOD GROWTH FUNCTIONS AND LINEAR YIELDS 

Models with linear yields were biologically motivated by (TJ [5l [6] who noticed the existence of 
limit cycles for some values of the parameters. The rigorous mathematical study was given in [18] . 
Consider the particular case of (0, where the growth functions qi(S) are given by (jTTJ), and the 
yields Vi(S) = Pi(S) / fi(S) are linear 

y i (S)=Y i (l + c i S) 

where Yi > and a > 0. System (|7J|, with D = 1 and S° = 1, takes the form 



(45) 



S" = 1 - S - 



a,S 



N 



iS 



^h+sYiii + asy 

2 — 1 ' 



h + S 



D, 



■N. 



The break-even concentrations Xi are given by (|3ip . In this section we give analytical conditions 
on the parameters of (|45[) so that conditions (fl"4")) . (TT5|) and (TT5|) are satisfied and Theorem II .21 can 
be applied. We need the following technical result. 




= 27bc 





Figure 3. On the left: the definition of the function c cr n(b). For each b < 1, 
the functions y = [c(l — 6) — 1] (in green) and y = 276c 2 (in red) intersect for 
c = c cr it(b). On the right, the numerical plots (in red) of the function c = c cr it(b), 
showing the quick increasing of c cr it (b) with b. 
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Lemma 6.1. The function Q(S) = ~ X + X + c ) j s decreasing over [0, 1] if and only if 

[c(l - 6) - if < 276c 2 . 

This condition is equivalent to either b > 1 or b < 1 and c < c cr it(b), where c cr it(b) is the positive 
zero o/[c(l - b) - l] 3 = 276c 2 . 

Proof Since 

the function Q(S') has an inflexion point for 5 = (6/c)a. The function Q(S') is nonincreasing over 
[0, 1] if and only if its derivative at the inflexion point is nonpositive, that is, P' I (6/c)s ) < 0. 

Straightforward computations show that this condition is equivalent to [c(l — b) — 1] < 276c 2 . If 
6 > 1 then the first term of the inequality is negative and hence the inequality if satisfied for all 
c > 0. If 6 < 1, then the inequality is satisfied if and only if c < c cr i t (b), see Fig. [3] The expression 
of c cr it(b) can be obtained by Cardan formulas. Notice that c cr i t (0) = 1 and c cr u{b) is quickly 
increasing with 6. 

Theorem 6.2. Assume that 

(46) Ai < A 2 < ••• < X N , Ai < 1, 

(47) either b\>\ or for each i > 1 satisfying Ai < 1, Ci < c cr it(t>i). 

Then the equilibrium _E* is globally asymptotically stable for J^5[ ) with respect to the interior of the 
positive cone. 

Proof 'Let us prove that (fT4")) . (|T5|) and (|2"4")l hold. The Monod function fi(S) is increasing. Hence, 
(fT4l holds. The function P^S) is 

= (I-5)(6 1 + 5)(l + c 1 5) ^ 

By Lemma [67X1 it is decreasing if and only if either 6i > 1 or 6i < I and c\ < c cr u(bi). Hence, (fT5j) 
holds. For each i > 2, the function gi(5) defined by (|23|) is 

(q\ - M^j 1 ~ ^ - y» «j -a g- Aj 

/i(6)Pi(6) ajOi-AA-Ai 

where Qj(5) = (1 ~ S)(bl + S)(1+c ' S) . Assume that 07} holds. By Lemma O the function Qi(S) is 
decreasing. Therefore, 

min Qi(S) = Qi{X\) > Q%(\i) = max Qi(S). 

0<S<Ai Xi<S<l 

Since Ai < Ai, the function S i-» gZ\ i s increasing. Therefore, 

j — Ai Ai 1 — Ai ib — Ai 

mm — — = t— > 1 > 



o<s<Ai 5 - Ai Ai 1-Ai a,<s<iS'-Ai 

Thus, 

Yi a{ — A S — \i Yi ai ~ Di 
mm 5j(5) > — mm — — min Q,-(5j > —Qi(\i), 

0<S<Ai a 4 CL\ — A 0<S<Ai 5 — Ai 0<5<Ai CLi a\ - A. 
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and 



?o\ ^ *i ai ~ Di S — Xi , . Yi Qi — Di 

max g,(i) < — max — max Qi(b) < ttQ(^)- 

\i<S<i ai a\ — D\ \i<s<\ o — Ai A;<s<i <Zj a\ — D\ 



Hence (l24l) holds. The result follows from Theorem 11.2 





Figure 4. The graphical depiction of conditions (|15|) and (|24|) in the proof of 
Theorem 16.21 The parameter values are c\ = 4, b% = 0.1, ai = 1, 62 = 0.15, 02 = 1 
Di = 0.6 and D2 = 0.55. Hence Ai = 0.15 and A2 — 0.18. On the left, the function 
Pi(S) (in green) and its derivative (in red) showing that Pi is decreasing and so 
fTB]) is satisfied. On the right, the function g2(S) for C2 = 5 (in red), C2 = 30 (in 
green) and C2 = 80 (in cyan). The condition (|2~41 is satisfied for c^ — 5 < c cr ii (0.11) 
and C2 = 30 > c cr it(0.1). It is not satisfied for c 2 = 80. Here c cr ,t(0.1) ~ 6.5, see 

Fig. m 



Theorem 16.21 extends Corollary 14. 31 which corresponds to the case where the yields are constant. 
Indeed, (|46)) is the same as ([32]) and, for constant yields, Cj = 0, so that the conditions (|47|) in 
Theorem 16.21 are satisfied. Notice that (|47|) is a sufficient and not necessary condition for the 
existence of a gap between the minimum of gi(S) over (0, Ai) and its maximum over [Ai, 1]. For 
instance, for the parameter values given in the caption of Fig. [H if C2 = 30 > 6.5 ~ c cr ,t(0.1), there 
exists such a gap, see Fig. SJ Therefore, Theorem 11.21 applies and predict that the equilibrium 
is GAS, even if Theorem 16.21 does not apply, since b\ = 0.1 < 1 and C2 = 30 > 6.5 ~ c cr ,t(0.1). 
However, for C2 = 80 > 6.5 ~ c cr ,t(0.1), there is no gap, see Fig. |4] Therefore, neither Theorem II .21 
nor Theorem 16.21 can be used. This model shows that in each particular example it is very easy to 
depict graphically the conditions (fT4")) , (|15l) and (fT6"|) of Theorem 11.21 and to see if Theorem 16.21 can 
be applied. 



7. Discussion 

We briefly survey some CEP results for (Ql. In the Monod case [IT] when the growth functions 
are of form (|lip . and assuming equal removal rates for S and all species, i.e. Di = D for i = 1 • • • N, 
Hsu, Hubbell and Waltman [11] proved the following CEP: every solution of ((4]) with positive initial 
condition satisfies 

lim S(t) = Ai, lim Xi(t) = Yi(5° - Ai), lim Xi (t) = 0, i > 2, 

under the additional assumption < Ai < S° and Ai < Ai for i — 2 • ■ • N. The predictions in |11| 
were tested in the laboratory by the experiments of Hansen and Hubbell [8J . Similar experiments 
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could be performed to test the predictions using microorganisms known to have variable yields. See 
[I] [U O [6] [18] for examples of such microorganisms. 

Hsu [9] used the Lyapunov function (|25[) to give a simple and elegant proof of the result in [11] 
for the case of different removal rates Di (see Corollary I4.3|) . Wolkowicz and Lu [26J , used the 
Lyapunov function (|17|) and extended the results of [9] by allowing more general growth functions. 
They identified a large class of growth functions, where the constant c, in (jTTJ) can always be found. 
Following Smith and Waltman [23J, I think that, despite the fact the Cj cannot be found for all 
growth functions, the work of Wolkowicz and Lu |26| represents a major step in the extension of 
the result of Hsu [5] to general growth functions. In the constant yield case, the CEP has been also 
proved under a variety of hypotheses by Armstrong and McGehee |3], Butler and Wolkowicz @], 
Wolkowicz and Xia [27] and Li [TJ. The hypotheses used in the papers 0H M LLIl LU HS1 W\ are 
summarized in Table 1 of [T3] . Lyapunov techniques in the chemostat were also used in [151 \T§\ . 

The variable yield case was considered, for n = 1,2 by Pilyugin and Waltman [18], with a 
particular interest to linear and quadratic yields, and by Huang, Zhu and Chang [12] . The general 
model |[?J) for N species, was considered by Arino, Pilyugin and Wolkowicz [2]. As noticed by 
these authors (see [2J, Section 3), in the case of constant yields (QJ, including the yield terms Yi 
in the substrate equation, as in Q, is mathematically equivalent to including the reciprocal in the 
microorganism equation instead. Indeed, (|4]) can be written 

JV 



s' = D[s a -s]-J2MS)^, 



x\ = [Yipi(S) - Di] Xi , z = 1 • - • JV, 

where Pi(S) = q y . . Since Yi are constant, the uptake terms Pi and growth terms q^ have the same 
monotonicity properties. Formally, the model ([7J5]) with variable yields yi(S) can be written 

(48) S ' = ^-^-Ys^S)^ 

x l = [qi(S) - Di]xi, i = 1 • • • N, 

where qi(S) are the growth functions, or equivalently, 

N 

(49) S' - D[S°-S\-^2pi(S)xi, 

x'i = [Vi(S)pi{S) -Di]xi, i = l---N, 

where Pi(S) = , s \ are the uptake functions. One of the important differences in the case that 
the yields are not constant is that the variable yield terms can lead to uptake and growth functions 
that have now different monotonicity properties. Moreover, in the case of constant yields, the yields 
terms Yi can be eliminated in (j4]) simply by passing to the variables m = YiXi- We obtain 

N 

S' = D[S°-S\- y £ j q i {S)u h 

8=1 

u'i = [qi(S) - Di]ui, i = l---N. 
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This change of variables means that we have changed the units in which the microorganisms were 
evaluated. There is no such trick to eliminate the yields terms in (|48|) of (f49|) . Therefore, care- 
ful attention to the interpretation of the yield terms resulting in their correct placements in the 
equations is necessary. For details and complements, the reader is referred to [2], Section 3. 

In the variable yield case, the CEP has been proved for ([7]), under some technical conditions 
on the function pi and qi . by Sari [2T] and Sari and Mazenc [22] (see Corollary 13.11 in Section 
[3] and Corollary 14.41 in Section 3]). It was also shown in [22] how Corollary 14.41 can be fruitfully 
used to analyze the stability properties of systems whose yield functions depend on the variable 

5. For instance, the CEP holds (see Corollary 5 in [52]) f° r the Monod model with constant 
yields replaced by either linear or quadratic functions of S, and under certain additional technical 
assumptions. Another application is given by the model of Pilyugin and Waltman |18| which was 
used to demonstrate that a periodic orbit was possible in the case of variable yield model. In this 
model, with two species, where one yield is constant and the other is cubic in S, it is shown in |22) 
that for some values of the parameters the CEP holds (see Corollary 6 in |22J). The problem of 
the existence of limit cycles in chemostat equations is not always well understood [20]. In the case 
of constant yields, numerical simulations of model (j3| have only displayed competitive exclusion. 
Our results concern also the case of variable yields, for which it is known [2] [T2J [TB] that more 
exotic dynamical behaviors, including limit cycles and chaos, are possible. Thus in the case of 
variable yields, it is of great importance to have criteria ensuring the global convergence to an 
equilibrium with at most one surviving species. The reader interested in biological motivations for 
the dependence of the yields on the substrate, may consult [2J [TB] and the references therein. 

The exclusion of periodic orbits in system © was obtained by Fiedler and Hsu (see Theorem 
1.1 in |7J) under the following conditions: for all 1 < i ^ j < N, and < S < 1 

(50) (S - Xi)fi(S) > 0, forS^X u 

(51) f^S) < 1 + fj{S) + (1 - S)^(S)/ ft -(S). 

Even if the result in [7J does not show the convergence to an equilibrium, it is interesting to compare 
the constraints on the functions fa and pi of [7J with our constraints. Actually, ([50)1 is stronger than 
(|14p . since our assumption requires (|50|) only for /i and allows the fi, for i ^ 1, to have other zeros 
than Xi in ]0, 1[. Let us compare the constraints on fi and pi imposed by the inequalities (|51|) to the 
constraints imposed by hypothesis ([TBI : notice that ([5"Tj) is a set of — ~ - conditions, while (fTB"]) 
is a set of at most N — 1 conditions. Moreover, the constants a, in the conditions ([TBI give more 
flexibility to these conditions. For instance, (|T6| are satisfied by arbitrary Monod growth functions 
and also by a large class of growth functions as it was shown in [22. 26]. On the other hand, (|5Tj) are 
not satisfied by arbitrary Monod functions, see formulas (6.10) and (6.11) in [7j. Hence the result 
in [7J does not recover the CEP, even in the classical and well established case of Monod functions 
and equal removal rates [3] . However our theorem recovers a lot of results of the existing literature. 
The beautiful geometrical techniques in [7J, whose purpose is to extend Bendixon-Dulac criterion 
to higher dimension are innovative and it seems likely that [7j is destined to be cited more often for 
these innovative techniques rather than for the specific results that are proved for the CEP. 

For the purpose of comparison between our result and the result of Fiedler and Hsu [JJ, we just 
mentioned two caveats on Theorem 1.1 in [7J: first, this theorem does not recover many of the 
biologically interesting classical examples where the CEP is known to hold, and second, it does 
not prove the convergence to an equilibrium. These caveats were already mentioned in [Jj , Section 

6. Another caveat must be signaled. Fiedler and Hsu claimed (see [7J, Section 6) that, in the 
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case N = 1 of a single species, condition (|5Tj) holds trivially and there is no periodic orbit for 
system (f37j) . It should be noticed that condition ([50)) is not sufficient to exclude periodic orbits. Of 
course, if p'(S) > 0, then the Bendixon-Dulac criterion can be applied to exclude periodic orbits 
(see Section [5]). This assumption on the monotonicity of p is not explicitly stated in Theorem 1.1 in 
[7j. Moreover, the condition p'(S) > would not be satisfactory from the biological point of view. 
Indeed, a variable yield term y(S) = q(S)/p(S) can lead to nonmonotone uptake term p(S) even if 
the growth term q(S) is monotone (see Section [SJ. 

Fiedler and Hsu, see Section 6 in [7J, claimed that the construction of Lyapunov functions in 
[21 HH [251 [27] strictly depends on the proportionality © required in equations Q . We show how 
the Lyapunov function used by Hsu himself [9] for the Monod case, more than thirty years ago, 
can be extended to the case of (TT]), where growth rates are not required to be proportional to food 
uptake (see Theorem 14. 11 in Section 0|. For that reason, the direct proof of Theorem 14. 1[ using the 
extension of the Lyapunov function of Hsu |9J, seems to be interesting in itself. Thus, I decided to 
give Theorem 14.11 and its direct proof, despite the fact that this theorem is a corollary of Theorem 
L2j(see Proposition 14.21 in Section SJ). 

We list some references to the existing literature which inspired our approach. The Lyapunov 
function (I2oT) used in the proof of Theorem 14.11 was introduced in [22J as an extension of the 
Lyapunov function (|25[) that Hsu used in [9] in the Monod case (Theorem 3.3 in [9]). In the case 
of one species, this Lyapunov function is equal to the function (|4"2j) used by Pilyugin and Waltman 
(Lemma 2.3 in |18]L as shown in Section [5j It is also a multiple of the Lyapunov function that 
Ballyk, Lu, Wolkowicz and Xia used in [27], page 1039 or [25], Section 3.3 (see Section 3.2 in 
[22]). The Lyapunov function (fT5)) used in the proof of Theorem 11.21 was introduced in [2T] as an 
extension of the Lyapunov function (TT71) that Wolkowicz and Lu used in [26 in the constant yields 
case (Theorem 2.3 in [26]). In the case of one species, it is a multiple of the Lyapunov function (|40|) 
used by Arino, Pilyugin and Wolkowicz (Theorem 2.11 in [2]), as shown in Section [5] 

In this work, we have analyzed a general model of the chemostat with several species competing 
for a substrate, under the assumption that uptake rates and growth rates are not proportional. 
Each species is characterized by its specific growth rate, its specific removal rate, and its variable 
yield. Our study reveals that the CEP holds for a large class of systems: the species with the 
smallest break-even concentration can be the winner of the competition if some supplementary 
conditions, involving the uptake and growth functions are satisfied. Hence, even if the break-even 
concentration are depending only on the growth rates and not on the yields functions, the issue 
of competition depends really on the yield functions. For instance, if one on the species exhibits 
a linear yield, and if the parameter in the yield is enlarged, then the equilibrium, where only the 
winning species survives, can be destabilized, and oscillatory coexistence of more than one species 
becomes possible. 
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